#include "LinearSystem.h"

void solveLinearSystem(compressed_vector& x, compressed_matrix& A, compressed_vector& b)
{
	boost::numeric::ublas::permutation_matrix<double> P(b.size());
	boost::numeric::ublas::lu_factorize(A, P);
	x = b;
	boost::numeric::ublas::lu_substitute(A, P, x);
}
